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We study planar two-dimensional quantum systems on a lattice whose Hamiltonian is a sum 
of local commuting projectors of bounded range. We consider whether or not such a system has 
a zero energy ground state. To do this, we consider the problem as a one-dimensional problem, 
grouping all sites along a column into "supersites" ; using C*-algebraic methods^, we can solve this 
problem if we can characterize the central elements of the interaction algebra on these supersite. 
Unfortunately, these central elements may be very complex, making brute force impractical. Instead, 
we show a characterization of these elements in terms of matrix product operators with bounded 
bond dimension. This bound can be interpreted as a bound on the number of particle types in lattice 
theories with bounded Hilbert space dimension on each site. Topological order in this approach is 
related to the existence of certain central elements which cannot be "broken" into smaller pieces 
without creating an end excitation. Using this bound on bond dimension, we prove that several 
special cases of this problem are in NP, and we give part of a proof that the general case is in NP. 
Further, we characterize central elements that appear in certain specific models, including toric code 
and Levin- Wen models, as either product operators in the Abelian case or matrix product operators 
with low bond dimension in the non- Abelian case; this matrix product operator representation may 
have practical application in engineering the complicated multi-spin interactions in the Levin- Wen 
models. 



The subject of Hamiltonian complexity theory is de- 
voted to the study of the computational complexity of 
various problems in quantum many-body physics. A gen- 
eral framework for this problem is as follows. We con- 
sider a quantum system whose Hilbert space is the ten- 
sor product of N different Hilbert spaces. Each of these 
N Hilbert spaces is referred to as the Hilbert space of 
a "site". We are interested in the case that each site 
has Hilbert space dimension that is poly(N) (indeed, in 
many practical settings it is 0(1)). The Hamiltonian 
will be a sum of at most poly(N) terms, where each term 
in the Hamiltonian acting on at most O(l) sites. Of- 
ten, a locality condition is imposed on these terms in the 
Hamiltonian: there is some given graph and each term 
in the Hamiltonian only acts on a set of sites which has 
small diameter with respect to the graph metric (follow- 
ing physics terminology, the term "site" refers both to a 
vertex of a graph as well as to a Hilbert space associated 
to that vertex, and also following physics terminology 
we refer to edges as "bonds"). Hamiltonian complexity 
theory then addresses the complexity of answering vari- 
ous questions about this Hamiltonian; in this paper we 
address questions about ground state energy, but other 
problems including time dynamics and correlations have 
been considered by others. 

One basic result is that determining whether or not a 
quantum many-body system has a ground state energy 
less than or equal to some given value is complete for the 
complexity class QMA-i; this contrasts with the classical 
case where the analogous decision problem for the ground 
state energy of a classical system is complete for NP. One 
important area of research relates to the different promise 
gaps in these settings: the decision problem in the quan- 
tum context is phrased as deciding whether the energy 



is less than some value Eq or greater than some value 
Ei where the difference in energies is called the "promise 
gap" . Roughly speaking, one can think of this as requir- 
ing one to approximate the ground state energy to an 
accuracy given by the promise gap. Thus far, there is a 
big difference in the promise gaps in the quantum and 
classical settings. The quantum result is known to hold 
only for promise gaps which arc l/poly(N) (though with 
these polynomially small gaps the result has been ex- 
tended down even to one-dimension^) , while the classical 
result is known to hold for promise gaps which are of or- 
der the norm of the Hamiltonian itself as follows from the 
PCP theorem^. While there has been much effort in try- 
ing to prove or disprove the quantum PCP conjecture^—, 
thus far the problem is open as to how large the promise 
gap can be made in the quantum setting to still have a 
QMA-complete problem. 

Another important area of interest is the study of 
quantum Hamiltonians which are a sum of commuting 
projectors. This set of quantum Hamiltonians lies in be- 
tween the classical case and the quantum case. Thus far 
all QMA-completeness results for quantum Hamiltonians 
construct Hamiltonians which are not in this set of com- 
muting projector Hamiltonians. One nice feature of a 
commuting projector Hamiltonian is that all the eigen- 
values arc integers, so that the issue of a promise gap docs 
not arise: that is, given that the Hamiltonian is a sum 
of commuting projectors one knows that if the energy is 
greater than some given Eq then it is at least [£"0] + 1- 

Now consider the difference between the commuting 
projector case and the classical case. In the classical 
case, not only do the individual terms in the Hamiltonian 
commute, they are all diagonal in a product basis and the 
ground state of the Hamiltonian is a product state. In 
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the case of commuting projectors, this may not be possi- 
ble. Indeed, there exist Hamiltonians whose ground state 
cannot be turned into a product state by a local quan- 
tum circuit with depth and range small compared to the 
system size. Prototypical examples of such Hamiltonians 
in two dimensions include the toric code^ and Levin- Wen 
models^. Such ground states are called topologically or- 
dered states, and this definition that a state is topologi- 
cally ordered if it cannot be turned into a product state 
by a local quantum circuit is a circuit definition of topo- 
logical order (see^ among others). On the other hand, 
other commuting projector Hamiltonians do have topo- 
logically trivial ground states, meaning that the ground 
state can be turned into a product state by a quantum 
circuit of small depth and range. One example is when 
the Hamiltonian is a sum of terms acting only on two 
sites at a time^ (see also^). 

Such topologically trivial states can act as a classical 
witness for the existence of a low energy eigenstate of the 
Hamiltonian. Given a state ?p which is given by a circuit 
of depth O(l) acting on some given product state, with 
each unitary in the circuit acting on at most O(l) sites, 
one can compute the reduced density matrix of ip on any 
set of O(l) sites in a time that is poly(N); hence, one 
can compute the expectation value of each term in the 
Hamiltonian in a time that is poly(N) and so one can 
compute the energy of the state ip in a time poly(N). So, 
if such a Hamiltonian has a topologically trivial ground 
state, there is a classical witness which can be check in 
polynomial time. 

In this paper we show that for several special cases of 
commuting projector Hamiltonians in two dimensions as 
defined below, the problem of determining whether there 
is a zero energy ground state is in NP (the special cases 
are defined below; we only sketch part of a proof for the 
general case, leaving the rest for a future work). Note 
that this includes Hamiltonians whose ground state is 
topologically ordered. Hence, our classical witness can- 
not simply be a description of a trivial state using a quan- 
tum circuit. Instead, our classical witness will include 
certain central elements of certain interaction algebras 
as discussed below. These central elements may include 
"string-like operators" , as described below, characteristic 
of systems with topological order. Thus the key step in 
this result is characterizing possible string-like operators 
which may appear. These operators arc called "string- 
like" because they may have support on a column of sites 
in the two-dimensional lattice (the number of such sites 
will be of order L, for a two-dimensional system of size L- 
by-L). To specify an arbitrary operator with support on 
that many sites requires specifying exponentially many 
(in L) different matrix elements, and hence even if we 
only specify each matrix element by approximating it 
with polynomially many bits, the description of such an 
arbitrary operator would require exponential resources. 
Hence, we will need to give a better characterization of 
these operators. 

Note that it is known that the analogous problem for 



classical Hamiltonian is NP-completeii and hence this 
quantum problem is also NP-complete as it contains the 
classical problem. Thus, it should come as no surprise 
that we will not show how to find the witness efficiently. 

We begin by defining the problem of two dimensional 
commuting local Hamiltonian, which we call 2DCLH, in 
section U and show how various different systems can be 
cast in the form of our problem definition. Then, in sec- 
tion [H] we review the C*-algebraic approach to commut- 
ing Hamiltonians and the role of central elements in the 
interaction algebra. We also define the ideas of "propa- 
gation" , "breaking" , and "masking" which are used later. 
We show how complicated central elements can arise in 
even classical models, but for such models we show how 
these elements can be simplified by a "refinement" to a 
problem in a which the central element is a product of 
operators on a single sites, each such operator commuting 
with the Hamiltonian. In scction lHTl we discuss an exam- 
ple of central elements that appear in a simple toric code 
model with topological order. We show how one cannot 
simplify the central elements in this case to a product of 
operators on single sites which commute with the Hamil- 
tonian. Indeed one cannot even "break the string in two" , 
as such a breaking would create endpoint excitations; this 
is defined later and is shown using central elements of a 
different interaction algebra which measure topological 
charge. In section IIV1 we consider a restricted class of 
Hamiltonian where no such elements measuring topolog- 
ical charge do not exist, and show how conversely one 
can "break" the string-like operators when these terms 
are absent and in section [Vj we show that such Hamil- 
tonians have trivial ground states and hence there exists 
a witness to the existence of a zero energy ground state 
for these Hamiltonians. Finally, in section IVI1 we con- 
sider the problem of constructing a witness to the zero 
energy ground state for arbitrary instances of 2DCLH. 
We prove that such a witness exists for another special 
case, and partially sketch the general case; the full proof 
of the general case will appear elsewhere. 



I. PROBLEM DEFINITION 

We now define the problem 2DCLH. We consider a 
square lattice, with size L-hy-L. We refer to one direc- 
tion of this L-hy-L lattice as the vertical direction and 
one direction as the horizontal direction. The sites will 
be described by horizontal and vertical coordinates each 
ranging from 0, L — 1. 

Let N = L 2 be the total number of sites. On each 
site we have a Hilbcrt space with dimension D at most 
polynomial in N. The Hamiltonian that was consider is 

h = J2(i-Pz) = Y,Qz, (i) 

z z 

where Qz = 1 — Pz and the Pz are projection operators 
supported on a set of sites Z . The sum is over sets Z 
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that consist of four sites around a plaquette. Wc have 
{Qy,Qz}=0 (2) 

for all Y, Z. 

Note: throughout this paper, if we write an operator 
Oa, where A is a set of sites, it is assumed that the 
support of Oa is the set A unless we specify otherwise. 
Similarly, if we write Oij where i,j are sites, then the sup- 
port is the set {i, j}, or if we write Oijk where k are 
sites then the support is the set {i, j, k}, and so on. This 
convention applies also in the case of operators with addi- 
tional superscripts or subscripts: for example, wc might 
define several operators 0\, O a , O a , ... all supported on 
the set of sites A. When we do this, we will usually in- 
troduce an index a taking a value from some finite set 
and write a given operator from that last 0\ , 0\ , ... as 
u A . 

The problem 2DCLH is to determine whether or not 
H has a zero energy ground state. Throughout this pa- 
per, we gloss over details of real number arithmetic, but 
we make a few comments about this now. Suppose, first, 
that one wants to check that a Hamiltonian H really 
is a sum of commuting projectors. An operator Qz is 
described by several real numbers, namely the matrix el- 
ements of the operator. Using a polynomial number of 
bits to approximate each of these numbers, one can check 
that two projectors commute up to exponentially small 
error and one can check that each operator is a projector 
up to exponentially small error. While we will describe 
our result later as showing that there is a classical wit- 
ness for a state with zero energy given that the Pz are 
projectors and commute exactly, the witness will include 
some real numbers. If we use a polynomial number of 
bits to specify these real numbers, what the verifier will 
be able to verify is that there is a state with energy less 
than some quantity which is smaller than an exponential 
of a polynomial in N (rather than verifying the existence 
of a state with strictly zero energy). We will omit the 
details of this, though, and instead describe a witness 
specified by some exact real numbers. However, given a 
promise that the Hamiltonian actually is a sum of com- 
muting projectors, if one proves that there is a state with 
energy less than some exponentially small quantity, the 
existence of a zero energy state follows immediately, since 
the eigenvalues of the Hamiltonian must then all be non- 
negative integers. 

Note that the restriction to projectors is not a serious 
restriction. Suppose we have a Hamiltonian H — hz 
where [hy , hz] = for all Y, Z but the hz are not nec- 
essarily projectors. Then, the ground state of H can 
be chosen as an eigenvector of each hz with eigenvalue 
Xz with the ground state of H having energy ^2 z Xz- 
Then, we can define a Hamiltonian which is a sum over 
commuting projectors by defining a new Hamiltonian 
H' = Y^zO-~ Pz) where Pz projects onto the eigenspace 
of hz with eigenvalue Xz- Then H' has a zero energy 
ground state if and only if there is a state which indeed 
is an eigenvector of every hz with eigenvalue Xz- So, the 



restriction to terms which are projectors in the Hamilto- 
nian is not an important restriction. 

While some other features of our problem may also 
seem like arbitrary choices (for example, a square lattice 
rather than hexagonal or triangular or other lattice, or 
the restriction to a planar problem rather than a problem 
on a cylinder or other topology), the problem definition 
is actually not that restrictive. Consider any other two 
dimensional lattice with bounded range interactions on 
that lattice. We can group several of the sites in the orig- 
inal lattice "supersites" such that the resulting lattice of 
supersites is a square lattice and such that each interac- 
tion term in the original Hamiltonian acts only on sites 
in the same plaquette. Adding all the interaction terms 
in a given plaquette together, we get a sum of commuting 
terms H — ^2 Z hz as in the above paragraph. Note that 
this grouping into supersites will entail some overhead, 
in that the dimension of each supersite will be larger 
than the original sites, but if the original interactions are 
bounded range then this will lead to only an O(l) factor 
increase. 

Also, the restriction to planar two dimensional sys- 
tems is not so serious. In subsection IIII Bl wc show how 
to cast Hamtilonians in cylindrical or spherical geometry 
as problems with planar interaction, by "squashing the 
geometry flat" . This leads to some overhead, in that the 
dimension D on each site of the resulting planar problem 
is larger than that of the original problem, but for cylin- 
der and sphere (and other geometries with genus 0(1)) 
this gives only an increase in the dimension to D°^\ 



II. C*- ALGEBRAIC METHODS AND 
SUPERSITES 

We begin with a review of the method for finding 
trivial ground states of commuting projector Hamiltoni- 
ans which are a sum of two-site and one-site interaction 
terms. This is based on results o£i£ and see also^i^. We 
then discuss a grouping into supersites that we will use 
for the 2DCLH problem. 



A. C* Methods 

A key role is played by the concept of an interaction 
algebra, which originates in Ref. [l3l 

Definition 1. Given an operator O and a set X we de- 
fine the interaction algebra of O on X to be the algebra 
supported on X generated by all operators of the form 
tr-y (OQ-y ) where the trace is over all sites in the com- 
plement of X and Q-^ is any operator supported on the 
complement of X . 

More generally we can also define the interaction al- 
gebra of an algebra: given an algebra A, we define the 
interaction algebra of A on X to be the algebra supported 
on X generated by all operators of the form try(OQx) 
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where the trace is over all sites in the complement of X 
and Qjr is any operator supported on the complement of 
X and where O is any operator in A. 

Finally, whenever a Hilbert space 7-L decomposes as 
% = T~i\®H.2, we can define the interaction algebra of an 
operator O on a space Hi to be the algebra generated by 
all operators tr2(0Q2), where the trace is over space H.2 
and Q2 is any operator supported on H.2- We can define 
the interaction algebra of an algebra on Hi similarly. 

We also define 

Definition 2. Given an algebra, we say that the projec- 
tors P a are the minimal projectors which generate the 
center of the algebra if P a are projectors for each a that 
generate the center of the algebra so that P a pP = for 

We write a two-body Hamiltonian as 

// E fl 'J+E%. ( 3 ) 

<i,j> i 

where Hij acts only on sites i,j and H^, acts only on 
site i. The Hij and H^i need not be projectors. The 
notation Yl^j 9> indicates a summation over all pairs of 
sites i and j that are neighbors in some given graph; in 
the next subsection, we will talk about a grouping that 
reduces certain problems to the case where the graph is 
a line, while in this subsection the graph is not specified. 

Suppose all the various terms Hij and Hi commute 
with each other. Consider any site i. Let A 1 ^ be the 
interaction algebra of Hij on i. The interaction alge- 
bras A l i , A lk commute for j ^ k and also these algebras 
commute with iZj,j. 

Let Hi denote the Hilbert space on site i. We assume 
that Hi is finite dimensional. Then, it is a fact from C*- 
algcbra that we can decompose Hi into a direct sum of 
Hilbert spaces , 

h, = 0h; (, » (4) 

a(i) 

where a(i) is some discrete index, and then further de- 
compose each such Hilbert space into a tensor prod- 
uct of spaces %°^] (where the product ranges over j that 
neighbor i) tensor producted with space "H"\ so that 

Hi = Hf ] = {u«f ® (g) WSS) , (5) 

a(i) a(i) <j,i> 

where the product is over j that neighbor i, such that 
each operator Hij can be decomposed as 

H ^ = E P^Pf^lf ' aU \ (6) 

a(i),a(j) 

where is the operator on Hi which projects onto 

and Hlf Ml) acts on the subspace of Uf* ®uf' ] 



given by ® Hj-ll an d such that Hij can be decom- 
posed as 

a(i) 

where H^f acts only on T-i^f 1 ■ 

The operators Pf* commute with each other and com- 
mute with the Hamiltonian for any i and any a. The P" 
can be taken to be the minimal projectors which gener- 
ate the center of the algebra generated by the interaction 
algebras on i of J^j f° r a U 3- Thus, we can there is 
a ground state which is an eigenstatc of all of these op- 
erators. We find it useful to define an "effective classical 
Hamiltonian" : 

Definition 3. Given a Hamiltonian which is a sum of 
commuting terms, each acting on at most two sites, we 
define the effective classical Hamiltonian H e '^ as 
follows. We will give two definitions, one as a function 
from a discrete set to the real numbers, while the second 
definition will be as an operator. One can straightfor- 
wardly translate between these definitions and it should 
be clear in context which definition we mean. 

Let P a (i) be the projectors onto central elements de- 
fined above. Let 

H eff (a(l),a(2), ...,a(N)) = min^^^^paH^^, Hip). 

(8) 

That is, H e ff the minimum over all states tp, with norm 
|V'| = 1 and such that P a ^i' = ip of (ip, Hip). 
Similarly, we can define a quantum operator 

H eff (9) 

N 

£ (j[P^)H^(a(l),a(2),...,a(N)). 

a(l),a(2),...,a(N) i=l 

We noted that there is a ground state which is an eigen- 
vector of all the operators P°, so there is a ground state 
ip such that P^^ijj = ip for all i for some given choice 
of the indices a(i). Given the decomposition ([6]), we can 
construct such a state ^ which has simple entanglement 
properties: this state tp is a product of states ipij and 

states ipi, where ipi.j is in the space ® and^i 

is in Tift . As shown ir£, such a state can be generated 
by a quantum circuit of small range and depth acting on 
a product state. Note that given a Hamiltonian H , if we 
can determine a choice of a(l), a(N) which minimzes 
H e ff(a(l), a(N) then we can construct the states tpij 
in a time polynomial in TV (again, we gloss over the de- 
tails of the real number arithmetic required to do this) 
since each such state can be obtained by diagonalizing 

a Hamiltonian acting only on the space H"_£j Cg> 
which has dimension at most polynomial in N. 
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FIG. 1: Illustration of grouping. Smalled filled circles are 
sites. Vertical coordinates are shown ranging from to L — 1. 
3 of the columns are shown with horizontal coordinates 0, 1, 2. 
The figure below the arrow represents all sites in a single 
column grouped into a supersites with dimension D L . 

B. Grouping into Supersites 

The discussion above applies to Hamiltonians which 
are a sum of two-body terms. However, our Hamiltonian 
Eq. (JT|) is not a sum of two-body terms since each projec- 
tor may act on four sites in a plaquette. See Fig. [TJ We 
group all sites in a vertical column into a single super- 
site. We label the columns by numbers C . In an abuse 
of notation, we use C both as a number between and 
L — 1 labelling the horizontal coordinate of these sites in 
the given column and we also use C as a set (so that a 
site i is in the set C if and only if i has horizontal co- 
ordinate C). Then, the Hamiltonian Eq. ([T]) is a sum of 
terms each acting on a pair of supersites. Thus, we have 
succeeded in turning the 2DCLH problem into a problem 
of two-body interactions. However, this comes at a large 
cost: the dimension of the Hilbert space on each super- 
site is D L , so it is exponentially large. Writing down 
an arbitrary operator on a single supersite requires ex- 
ponential resources. The main work of this paper will be 
to show how to deal with this, by better characterizing 
the central elements of the interaction algebra that can 
appear in this problem. 

We need some definitions now. Recall that given an 
algebra A, the center of the algebra is the set of elements 
of the algebra that commute with every other element 
of the algebra. These elements are also called central 
elements. 

Definition 4. The projectors Qz generate a subalgebra 
of the algebra of all possible operators on the D N dimen- 
sional Hilbert space. We call this subalgebra the term 



algebra. 

Given a set S , we say that the term algebra of S is 
the algebra generated by all Qz with Z C S. Given sets 
S,B, we define the interaction algebra of S on B to 

be the interaction algebra of the term algebra of S on B. 

We define that 

Definition 5. A projector P is consistent with a 
Hamiltonian H written in the form of Eq. [Ity if there 
is a zero energy ground state of H that is in the range 
of P. Note that for a projector to be consistent with a 
Hamiltonian, that Hamiltonian must have at least one 
zero energy state. 

A projector P' is a refinement of P if P' < P. P' is 
a consistent refinement of P if P 1 is consistent with 
the Hamiltonian H and P' is a refinement of P. 

Finally, we define that 

Definition 6. An operator Ox supported on a set X 
which commutes with the term algebra is breakable into 
sets Y\, Y2, ... if we can write 

°* = E(n°n) ( 10 ) 

a a 

where the sum is over a in some finite set ( when we write 
an operator as a sum of product of operators with a sum 
of this form we will suppress the set from which a is 
chosen when writing the sum, as it should be clear from 
context) and for some choice of operators Y such that 
Y is supported on set Y a and O y commutes with the 
term algebra for all a, a. The operators O y need not 
commute with Yb ; if they do not commute, then we re- 
quire that there be some choice of the ordering in Eq. U0\) 
for which the equation holds and different orderings can 
be chosen for different a in the sum. 

The reason for this last definition is that we will be 
interested in projectors that arc central elements of the 
interaction algebra that are not breakable into sets of 
small diameter or which cannot be refined into projec- 
tors which are breakable on sets of small diameter. Let 
us consider first a case with breakable projectors to illus- 
trate what can happen when the projectors which appear 
as central elements of the interaction algebra are break- 
able. Consider as an extreme case a classical model in 
two dimensions with nearest neighbor interactions, where 
by classical we mean that every interaction term is diag- 
onal in some product basis. Let this product basis be 
the product of the basis on each site with states labeled 
|1), |2), ...\D). Then, the interaction algebra of hc-i,c 
on column C includes central elements which are also di- 
agonal in this product basis. However, assuming that H 
has a zero energy ground state then in this case all of 
these central elements which are consistent with H have 
a consistent refinement to a projector P which is a prod- 
uct of projectors on single sites in the column with each 
such single site projector having rank 1, so that 

p c = n p* ( n ) 

iec 
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To see this, simply take any zero energy ground state of 
H which is a product state (such a state exists because 
the Hamiltonian is classical) and let Pi be the reduced 
density matrix of that state on site i. Each of these 
projectors Pc commutes with the term algebra. When 
the central elements are breakable in this fashion, we can 
use the projectors on small sets as a classical witness 
for the existence of a zero energy state. However, we 
will see in the next section that in certain systems with 
topological order there exist central elements which are 
not breakable into small sets as discussed in subsection 

MM 

C. Propagation 

Finally, we describe an iterative method of determin- 
ing whether a given instance of 2DCLH has a zero energy 
ground state. In order to verify that there is a zero en- 
ergy ground state, we wish to verify that tr(J| z Pz) > 1 
as this trace counts the number of zero energy states. 
We can compute this trace using the following iterative 
procedure. Define 

Pc.c+i= n p ^ ( i2 ) 

zeCuc+i 

where the product is over sets Z containde in the union 
of column C with column C + 1. For use later, define 

Pc,c>= n 

ZeCUC+lUC+lU...UC" 

for any pair of columns C, C with C < C 
Define pc by 

Po = I (14) 

and for C > define 

pc+i = tr c (pcPc,c+i J , (15) 
where the trace is over all sites in column C. Then, 

tr(J[Pz) =tv L -i{pL-i). (16) 

z 

The operators pc can be written as matrix product 
operators^—. However, the bond dimension required to 
write a given operator pc may grow exponentially with 
C, so the above method may take exponential time. To 
show that the problem is in NP, later we will construct 
a classical witness which will help reduce the bond di- 
mension of the matrix product operators that we need to 
consider. Note that each pc is in the interaction algebra 
of C — 1 U C on column C (recall that in our terminology 
this interaction algebra is the interaction algebra on col- 
umn C of the algebra generated by the Q z supported on 
the union of column C — 1 with column C). Things can 



be simplified to some extent, since we only care about 
the center of this algebra; that is, let Pg be the minimal 
projectors which generate the center of the interaction 
algebra of the term algebra on column C . That is, the 
center of this interaction algebra is the algebra generated 
by the Pg and tv(PgP§) = 5 aJj . Then defining 

<fc=E P ctr(Oc), (17) 

a 

we have that 

pc+i = ti c (^cPc,c+iPz), (18) 

We say that pc is the propagation of pc from C to 
C + 1. The propagation of pc from column C to C for 

C > C is given by tTc,c+i,...,C'-i\PcPc,C'\ Note that 
this is equivalent to propagating pc from C to C+l, then 
propagating the result to C + 2, and so on, until reaching 
column C. (Note that a similar idea of propagation was 
used ir>i£ as a heuristic method of finding ground states of 
classical Hamiltonians using a heuristic method to reduce 
the bond dimension.) 

We now define "masking" : 

Definition 7. We say that an operator Oc supported on 
a column C masks a given operator Oc—i supported on 
column C — 1 if 

trc(Oc-iPc-i,cO c Pc,c+i) (19) 
= a;tr c (o c Pc,c+i) 

for some scalar x. 

We say that Oc is a mask if it masks all operators 
Oc-i supported on column C — 1 (note that the constant 
x may depend upon the given operator Oc-i chosen). 

Note that the constant c in Eq. (fT^) may equal zero. 

Note that every minimal central element of the inter- 
action algebra of C — 1 U C on column C is a mask. Also 
every minimal central element of the interaction algebra 
of C— lUCUC+1 or C, C + l on column C is a mask. 

Masking can be used to construct a witness for the 
existence of a low energy state as follows. 

Lemma 1. Suppose that for each column C , there exists 
some projector Pc which is a mask ( these projectors Pc 
should not be confused with the projectors Pz used to de- 
fine H ; it should be clear by the subscript which projector 
is meant). Suppose that Y\ c Pc is consistent with H and 
suppose that 

[P C ,Pc-i,c}=0 (20) 

for all C with 2 < C < L. Suppose finally that each such 
Pc is represented as a matrix product operator with bond 
dimension that is poly(N). Then, given such projectors 
written as matrix product operators, one can verify the 
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existence of a zero energy state in a time polynomial in 
N. 

Further, even if such projectors Pc are given only for 
some subset of the columns C , with the projectors being 
given for the columns in the sequence Ci,C2,...,C n for 
some n with 1 < C\ < C2 < ••• < C n < L, such that 
Cj_|_i — Ci < 0(1) for all i and such that L — C n < O(l) 
and C\ — 1 < 0(1), one can still verify the existence of a 
zero energy state. 

Proof. As mentioned before, we gloss over details of real 
number arithmetic, assuming all calculations are done to 
infinite precision. 

We only describe the case in which projectors are given 
for every column C in detail; the second case can be re- 
duced to the first case by defining a new problems where 
columns 1, ...,Ci are grouped into a single column (that 
is, all sites in those columns with a given vertical coor- 
dinate are grouped into a single site in the new problem 
with the same vertical coordinate), columns Ci + 1, C2 
are grouped into a second column, and so on. Since 
Cj+i — Ci < O(l), the Hilbert space dimension of the 
new problem is D°^\ 

The verifier first checks that Eq. (f2"0)) holds. Since the 
Pc are given as matrix product operators, this can be 
done efficiently 

Since Yic P c 1S consistent with H, 
tr (ric Pc Tiz p z) > 0. Conversely, given that 
tr (llc Tlz ^z) > 0, the existence of a zero en- 
ergy state follows. So, the verifier just needs to check 
that tr(n c Pc Hz p z) > 0- However, 

L-l 

tr(n^n^) = tr ([i p c n p ^ c+ i) (21) 

C Z C C=l 

= tv{P 1 Px t2 P 2 P 2 ^..P L -xP L -x tL P L ), 

where the last line follows by Eq. ([20| . 

For each column C with 1 < C < L — 2, the verifier 
computes the propagation of Pc from C to C + 1 , calling 
the result pc+i- The result is a matrix product operator 
with bond dimension poly(N), and this calculation can 
be done efficiently (note that the bond dimension of the 
resulting operator is at most a factor bigger than 

the bond dimension of Pc ■ The verifier than propagates 
pc+iPc+i to C + 2 and compares to pc+2- The verifier 
checks that the propagation of pc+iPc+i to C + 2 is 
indeed equal to a constant times pc+2, as in Eq. (|19[) . 
Call the constant for which Eq. ([19]) holds cc+i- The 
verifier also checks that the constant cc+i is positive. 
Since 

L-l 

t1iP1P1.2P2P23-PL-1PL-1.LPL) = J] c c tr(P t _ 1 P i _ 1 , i ) 

C=2 

(22) 

checking positivity of each cc and of the trace on the 
right-hand side of Eq. (|22|) is equivalent to checking pos- 
itivity of tr(n c Pc Uz p z)- □ 



III. CENTRAL ELEMENTS IN TORIC CODE 
AND OTHER MODELS 

We now describe a model with central elements that 
cannot be broken into certain smaller sets. We begin with 
a description of a toric code model on a square lattice 
and describe the central elements and effective classical 
Hamiltonian. We then show the unbreakability. Next, 
we consider some other geometries. Finally, we discuss 
the central elements that can appear in Levin- Wen mod- 
els among others, where the central elements are matrix 
product operators. 

Consider the toric code HamiltonianX: we have a 
square lattice with spin- 1/2 degrees of freedom on each 
site of the lattice (often the degrees of freedom are placed 
on the bonds in the literature but we place them on the 
sites here). The Hamiltonian is a sum of commuting 
terms. Color the plaquettes of the square lattice alter- 
nately as light or dark plaquettes, like a checkerboard. 
The Hamiltonian is 

h =- e IK- e n<< ( 23 ) 

Plight ieP PedarkieP 

where the first sum is over light plaquettes P and the 
product is over sites i in the plaquette and the second 
sum is over dark plaquettes and the operators af , a? are 
Pauli spin operators on site i. Note that here we choose to 
consider Hamiltonians which are not a sum of commuting 
projectors in order to keep our notation closer to the 
literature; recall the discussion of this in section |U 

We will discuss this Hamiltonian on various geometries, 
namely square, cylinder, torus, and sphere geometries. 
We begin with the square geometry, as in section |TJ In 
the square geometry case, we will add some additional 
terms to the Hamtilonian on the boundaries: on each 
pair of neighboring sites i, j on either the top or bottom 
edge of the system (these two sites will be connected by 
a horizontal edge), we add the term — erf trj if those two 
sites are in a light plaquette. Now, group all sites along a 
single vertical column C into one "supersite" as described 
above. Define O c to be the product of erf over all sites i 
in a given column C: 

O c = H*f. (24) 

iec 

The operator O c is a central element of the interaction 
algebra of C U C — 1 on C and also of the interaction 
algebra of CUC + 1 on C for all C. Further, the effective 
classical Hamiltonian 

H eff = J2 h c f c + i, (25) 
c=o 

involves these central elements, where 

frgo+i = -OhOh+v (26) 
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FIG. 2: Lattice for the toric code with L = 6. Light and 
dark plaquettes are indicated, with light plaquettes in the 
four corners of the system. Interactions are as described in 
the text. 



This square lattice geometry shows the feature we have 
mentioned: the existence of central elements in the in- 
teraction algebra which cannot be broken into certain 
smaller sets, as we will show in the next subsection. The 
physical reason for this unbreakability is that breaking 
the operator will necessarily lead to the creation of an 
anyon excitation. 

However, this square lattice geometry also lacks one in- 
teresting feature: there is no term in the effective classi- 
cal Hamiltonian on either the first or last supersitc which 
would constrain the central element on that super- 
site to have any particular value. That is, this Hamilto- 
nian has two degenerate ground states corresponding to 
different values of Oq. We can add this interesting fea- 
ture to the Hamiltonian by adding additional interaction 
terms on the left and right edge of the system. On each 
pair of neighboring sites i,j on either the left or right 
edge of the system (these two sites will be connected by 
a vertical edge), we add the term — erf cj if those two sites 
are in a light plaquette. Further, we choose the geometry 
so that the four corners of the system are all light pla- 
quettes (this means that we choose L to be even). Note 
that the corner sites now have two such added interac- 
tion terms, — erf o"J, one coupling a corner site i to a site 
j neighboring it by a vertical bond and one coupling it to 
a site j neighboring it by a horizontal bond. See Fig. [2] 

By adding these terms to the Hamiltonian, we have 
now the effective classical Hamiltonian 

L-l 

H eff = J2 h tc + i-°o-° X L-n (27) 
c=o 

which now has a unique ground state that minimizes ev- 
ery term in the Hamiltonian separately. We could also 
have considered a slightly different Hamiltonian, by pick- 
ing one of the vertical edges on the left side of the sample 
and changing the sign of interaction term so that we add 
+af(jj rather than — crfcr?. In this case, the effective 



classical Hamiltonian is 

L-l 

H eff = £ h e Jf c+1 + 01 - Ol_ 1: (28) 
c=o 

which has no ground states that minimize every term in 
the Hamiltonian. 



A. Unbreakability of Central Element 

Recall that we each site has a horizontal and a verti- 
cal coordinate, each coordinate ranging from 0, L — 1. 
A column C contains all sites with a given horizontal 
coordinate. 

Consider the toric code Hamiltonian on the square lat- 
tice with these added operators —o~fo~j on the boundary. 
We now show that 

Lemma 2. Suppose that 2 < C < L — 3. The central 
element Ofj is not breakable into any sets Yi,Yz, ... which 
have the properties that 

1 For all a the set Y a does not contain any sites with 
horizontal coordinate or 1 or horizontal coordi- 
nate L — 2, L — 1. 

2 At least one of the following two properties hold: 

[2a]Y a does not contain any sites with vertical 
coordinate 0. 

[2b] Y a does not contain any sites with vertical 
coordinate L — 1. 

Remark: requirement (2) above can be roughly stated 
as Y a can be close to the top edge or to the bottom edge, 
but not to both edges while requirement (1) can be roughly 
stated as Y a cannot be close to either the left or right edge. 

Proof. Suppose instead Oq is breakable. Let M be the 
set of sites which have horizontal coordinate or L — 1, 
or have vertical coordinate or L — 1, or both. That 
is, M is the set of sites on the boundary of the system. 
Define the operator T to be the product of af over all 
sites i in M. 

The operator T is in the algebra generated by the Qz, 
because the product of Qz for all light plaquettes is equal 
to T. Indeed, this operator T is a central element of the 
interaction algebra of the term algebra on M . Since each 
operator O y commutes with each Qz-, 

[O£ a ,T}=0. (29) 

Note further that \O x c ,T] = 0. 

Define T* to be the product of erf over all i in M which 
have vertical coordinate less than or equal to L/2 and 
let T b be the product of af over all j in M which have 
vertical coordinate greater than L/2, so that 

T = T f T b . (30) 
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Note that 

{O x c ,T t } = {O x c ,T b } = 0. (31) 

For each a, at least one of property 2a or 2b holds. 
Suppose 2b holds. Then, [Oy a , T b ] = since the opera- 
tors O y and T b have disjoint support. By Eq. (|29|) we 
have 

Oy a T l T b = T l T h 0^ a (32) 
= T t O Ya T b . 

Since T b is an invertible operator, we can multiply the 
right-hand side of the above equation by (T 6 ) -1 to obtain 

[O Ya ,T t ] = 0. (33) 

Thus, O y commutes with both T* and T b . We can sim- 
ilarly show that Oy commutes with both T* and T b if 
2a holds. 

So, the product over all a of O y commutes with T* and 
T b for each a. However, Oq anti-commutes with both T 4 
and T b . Therefore, Oq cannot be equal to FJ O y . 

□ 

The key role in this lemma was played by the existence 
of T which is also a central element in an interaction 
algebra. We will show later in section [V] that if certain 
similar central elements do not appear, then there are 
topologically trivial (under a circuit definition) ground 
states. 



B. Other Geometries 

Now consider this toric code Hamiltonian on a cylin- 
der, torus, and sphere geometries. We can take the toric 
code Hamiltonian on such geometries and embed it into a 
Hamiltonian on a finite square lattice with open bound- 
ary conditions with still a small Hilbert space dimension 
on each site. Consider the case of a cylinder, first: sup- 
pose the system has L y sites in the direction with periodic 
boundary conditions and L x sites in the direction with 
open boundary conditions. Label the site on the cylinder 
with coordinates (x, y) with < x < L x and < y < L y . 
Then, define a system on a square lattice with length 
L x in the x-direction and L y /2 in the y-direction with 
coordinates (x,y) with < y < L y /2. Let each site on 
the square lattice have a four-dimensional Hilbert space, 
with site (x, y) on the square lattice containing the de- 
grees of freedom of sites (x, y) and (x, L y — 1 — y) on the 
cylinder. 

This procedure can the thought of as "squashing the 
cylinder flat" : consider the cylinder as a two-dimensional 
surface embedded in three dimensions and project it onto 
a two-dimensional plane. This leads to only a small in- 
crease in the Hilbert space dimension on each site. The 
same kind of procedure can be applied to the torus or 
sphere case. 



For the case of the torus Hamiltonian, there is a four- 
fold degenerate ground state, with different ground states 
having different expectation values for these central ele- 
ments of the interaction algebra, and so for this Hamil- 
tonian it might seem like considering these central ele- 
ments is not important for determining whether or not 
there is a zero energy ground state. Consider, however, 
the sphere case. In this case, there is a unique zero en- 
ergy ground state. Suppose, however, that we make a 
change to the Hamiltonian on the sphere, changing the 
sign of the interaction term on one light plaquette. Then, 
there are no zero energy ground states. If we change the 
interaction sign on another light plaquette, even on a 
light plaquette located far away from the first plaquette 
with a changed sign, then there is again a zero energy 
ground state. This effect is due to central elements in 
the interaction algebra. Suppose we squash the sphere 
flat so that the north pole is at the left-most edge of our 
one-dimensional chain of supersites and the south pole 
is at the right-most edge of our chain of supersites. If 
all plaqucttcs have the same interaction sign except pos- 
sibly for those at the north and south poles, then the 
effective classical Hamiltonian is similar to Eq. (|2"ij|) in 
that it constrains values of the central elements to be the 
same between neighboring supersites, but there are ad- 
ditional terms at the left-most and right-most supersites 
constraining the central elements on those sites to have 
a given value depending upon the sign of the interaction 
term on the corresponding plaquette, just as in Eq. (|27[) . 



C. More General Models 

In the toric code model, the central elements O^ is a 
product of local operators. Such a product of local oper- 
ators is a matrix product operator with bond dimension 
1. Using the square lattice toric code with boundaries 
conditions as above, the central elements <fic of Eq. (TTTj) 
are the operators 1 + O^ which can be written as matrix 
product operators with bond dimension 2. If we consider 
several copies of the toric code, we can increase the bond 
dimension still further: consider a system with /c-qubits 
on each site describing k copies of the square lattice toric 

code. Then the operator <f>c is IIa=i(l + ^c( a ))> wnere 
Oq{o) is the operator O^ acting on the a-th copy. Such 
an operator cj>c has bond dimension 2 k . 

However, while this toy model of several copies of the 
toric code shows a case where we can increase the bond 
dimension needed to represent <f>c as a matrix product 
operator, the central element 4>c is a sum of a small num- 
ber (in this case, 2 k ) different central elements, each of 
which is a product operator (for example, one such sen- 
tral elements is Oq(o), for any given a. In contrast, if we 
consider more general Levin- Wen models, there may be 
no central elements that are product operators, and in- 
stead all the central elements may be matrix product op- 
erators (it seems likely, but we have not shown this, that 
a Levin- Wen model will have central elements that are 
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product operators if and only if it is an Abelian model) . 

Conside a Levin- Wen model^ on a hexagonal lattice, 
with the degrees of freedom on the bonds of the lattice. 
This paper gives the plaquette operators as twelve- 
spin operstors, acting on the six bonds around a given 
hexagon as well as the six external bonds connected to 
that hexagon (here, p labels a given hexagon in the lattice 
and s is a particle type). The matrix elements of this 
operator Bf> are written in terms of the F-symbols of 
the theory. There is a natural generalization from this 
operator B s v to an operator Bf acting around any non- 
self-intersecting closed loop I. Let the loop contain n 
bonds and label the bonds along the loop as 1, 2, 3, n 
around the loop in a clockwise direction, and the external 
bonds connected to the loop as ei, e2, e2, e n , with bond 
ei in between bonds n and bond 1 and in general bond 
e a in between bonds a — 1 and a. The operator Bf does 
not change the state of the bond e\, ...,e„. The matrix 
elements of £?; are given by 

(1',2',3'j ...;ei,e 2 ,e 3 , ...\Bf |1, 2, 3, e 1 , e 2 , e 3 , ...} 

R s,l',2',3',... 
— 1,1,1, 3,... 

= Fp^Fpl^Fp^,.... (34) 

In the above equation, the a term e a in an F-symbol 
means the state of bond e a , while a term a represents 
the initial state of bond a and a term a' represents the 
final state of bond a. The last line in the above equation 
can be written as 

n 

with the subtraction being done mod n so that 1 — 1 = n. 

If there are nt yP es particle types in the theory (plus also 
the identity particle), then the operator Bf can be writ- 
ten as a matrix product operator of low bond dimension 
as follows. Before describing how to write it as a matrix 
product operator, note that a matrix product operator is, 
by definition, an operator that acts on a one-dimensional 
system (a line). The spins that we consider do not quite 
form a one-dimensional system, since the external bonds 
"branch off' of the hexagon. To deal with this, we will 
choose to combine the bonds e a and a into a single de- 
gree of freedom so that we can remain in the formalism 
of matrix product operators. To define a matrix product 
operator, we need to also introduce auxiliary degrees of 
freedom; let a a be the auxiliary degree of freedom that 
couples the operator on e a and a to the operator on e a +\ 
and a + 1 . Then we will show how to write 

B i = XX°"" Q1 °2 2Q1 03 3Q2 ...) (36) 
{<*} 

= tr {Q} (0 1 2 3 ...), 

where the sum on the second line is over all possible 
assignments to the auxiliary degrees of freedom a, and 



where On"""" 1 is an operator acting on e a and a that 
depends upon the choice of the auxiliary degrees of free- 
dom a a and a a -i- The trace on the last line is over the 
auxiliary degrees of freedom a and the product 0\02 - - is 
a product of the matrices (viewing O a as being a matrix 
of operators). 

We let the auxiliary degrees of freedom a a be made up 
of two indices fi a ,P' a - Each index P a ,P' a will take one of 
n>t yP es + 1 possible values corresponding to the different 
particle types, including the identity particle. We define 

O a " " "~ 1 a ~ 1 by its matrix elements: 

{a'-eapt^-^'-^a-ea) (37) 

X X TjieaiPa-l)* a 

This may be more clearly written as 

(a';e a \Of^'\a;e a ) (38) 

_ x X 

where we have set [i = /3 a _i, v = (3 a to reduce the number 
of subscripts and superscripts. 

Eq. ([57]) gives the operators to write Bf as a ma- 
trix product operator with bond dimension (nt yP es + l) 2 - 
These operators Bf are central elements of the interac- 
tion algebra on I of the term algebra. For the doubled 
Fibonacci model, for example, Bf has bond dimension 4. 



IV. BREAK ABILITY 

In section [THAI we used the existence of an operator T 
which was a central element in an interaction algebra to 
show show the unbreakability of certain operators. We 
now show a result that is in some sense a converse of this: 
we show that for models with certain properties of the 
interaction algebras, all operators which commute with 
the term algebra are breakable. 

Let B be some set of sites which is the boundary of a 
set S of sites, where the boundary of S is the set of sites 
i € S such that i e Z for some plaquette Z <£_ S. For 
example, see Fig. [3J or Fig. [4] The sets / and E defined 
in the figure are the "interior" and "exterior" of B, as 
illustrated, with IC\B = EC\B = $; that is. / = S \ B 
and E = A \ S, where A is the set of sites in the whole 
lattice. Note that every term Qz is supported on S or 
E U B. Note that it is possible for a term Qz to be 
supported on B (for example, if S is the set of sites in a 
given plaquette, so that B = S) and hence be supported 
on both I U B and E Li B. Let A 1 be the interaction 
algebra of S on B and let A E be the interaction algebra 
on B of the algebra generated by the Qz for Z gt S. 

Let %b be the Hilbert space on B. The algebras A E 
and A 1 commute. So, if at least one of these algebras 
does not contain central elements, then we can decom- 
pose TLb = Hb^e <& T~Lb^,i so that A E acts only on 
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FIG. 3: The sites of the lattice are on the intersections of the 
solid horizontal and vertical lines. The heavier dashed line 
represents the set B. B is a closed line, containing the 12 
sites around around a 4-by-4 square. The 4 sites inside the 
square represent I and E consists of all sites not in / or B. 
The heavy solid line is A. The set X is A \ I, so X consists 
of all the sites in A above the line B and also the intersection 
of A and B. The set X has been indicated by drawing light 
diagonal lines around the sites in that set. 

Hb^e and similarly A 1 acts only on Hb^i- This de- 
composition will be the condition under which we show 
brcakability below. 

Before showing brcakability of certain operators, we 
give a lemma which has some physical interpretation in 
terms of the absence of "anyon excitations" as discussed 
later. 

Lemma 3. Define some set S with boundary B as above. 
Define another set A which is an arbitrary set of sites. 
Define 

X = A\I. (39) 

(We illustrate a possible choice of such sets in Fig. \^) 

Suppose Hb decomposes into Hb^e ® Hb^i so that 
A E acts only on Hb^e and similarly A 1 acts only on 
Hb->i- Let O be any operator supported on A such that 
O commutes with all Qz with Z D E ^ 0. Then, we can 
write 

= J2(°xuBO a BU(A \ x) ), (40) 

a 

where the operators 0^ uB are supported on X U B and 
the operators B[Jl ^ A ^ x ^ are supported on B U (A \ X) = 
B U (A H I) and where XuB commutes with the term 
algebra. 

Proof. The set A is a subset of the union of the disjoint 
sets X\B, B, and AC\I. The Hilbert space H B on B can 
be decomposed as a product Hb^e ®Hb^i- Thus, the 
Hilbert space on the support of O can be decomposed as 
a product of two Hilbert spaces H ex t <8> TLint , with 

H ext = U x \b ® Ub^e (41) 



and 

H m t=n B ^i®Usm- (42) 
Hence, we can decompose O as a sum 

= Y.°e X tO? n u (43) 

a 

where O ext acts on H ex t and Oi n t acts on Hint, and fi- 
nally where we pick that the operators tr {iPf n t)^ 0^ n ^\ = 

for a 7^ p and similarly tr (f0^ xt )^ Og Xt J — for a ^ (3. 

That is, we choose the operators to be orthogonal with 
respect to the Hilbert-Schmidt inner product. Such a 
choice is possible using a singular value decomposition^. 

We claim that for all a the support of 0^ xt is a subset 
of XUB and the support of Of nt is a subset of BU(Anl). 
To see this, note that O commutes with any operator sup- 
ported on E\X, and so using the orthogonality property 
of the Of nt , each operator 0" xt commutes with any op- 
erator supported on E \ X. Hence, 0" xt is supported on 
the X U B. The claim for Of nt follows similarly. 

Therefore, Eq. (|4"3"| gives us a decomposition of O in 
the form of Eq. (|4*U)) . We just need to show that 0° xt 
commutes with Qz for all Z . If Z fl I ^ 0, this follows 
immediately since Qz acts on Hint m that case. However, 
if Z PI / = 0, then Qz acts only on % ex t and not on Hint 
and so, given that [Qz,0] = and that [Qz,Of nt ] = 

and given the orthogonality property of the Of ntl it 
follows that [Qz, 0^ xt ] — for all a. 

This completes the proof. We note that the last para- 
graph is simply another application of the concept of an 
interaction algebra. Let the interaction algebra of O on 
H ex t be the algebra generated by the 0" xt . Since O com- 
mutes with Qz and Qz does not act on Hint, the inter- 
action algebra commutes with Qz- □ 

This last lemma [3] has a physical interpretation that 
there are no anyon degrees of freedom inside /. For a 
theory like the toric code which has anyons, there are 
operators which are string-like which create excitations 
at their endpoints. For example, consider the product 
of erf over all i in a column C such that the vertical 
coordinate of i is less than or equal to y for some y < L. 
This operator does not commute with the Hamiltonian, 
but it commutes with every term except for those terms 
Qz with Z acting on the endpoint of this string. Indeed, 
in the toric code it commutes with all but one term in the 
Hamiltonian, so that acting on the ground state it creates 
a state with energy 1. One can show that this operator 
erf cannot be decomposed in the form of Eq. (|4T)1) for 
any set B centered on the endpoint of the string with 
B small enough that B does not touch the boundary of 
the system; we omit the proof of this statement as it is 
very similar to the proof of lemma [5] and is based on 
considering commutation with a certain central element 
of an interaction algebra. 

Similarly, one can consider the product of erf over all 

1 in column C such that the vertical coordinate of i is 
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FIG. 4: The sites of the lattice are on the intersections of 
the solid horizontal and vertical lines. The heavier dashed 
line represents the set B. B is a closed line, containing the 
12 sites around around a 4-by-4 square. The 4 sites inside 
the square represent I and E consists of all sites not in I or 
B. The heavy solid line is A; one bond in this line has been 
draw even thicker, so that sites in A above the middle of this 
line are in A\ and sites in A below the middle of this line are 
in A-2- The sets X, Y have been indicated by drawing light 
diagonal lines around the sites in that set, with X consisting 
up the upper set of such points and Y consisting of the lower 
set of such points. 

between x,y for some 1 < x < y < L. This operator 
creates a pair of excitations at the ends of the string. 
It commutes with all but two terms in the Hamilto- 
nian, one near one end of the string and one near the 
other end. One can similarly show that this operator 
cannot be decomposed into a sum of products of oper- 
ators, J2 a 0"0"0", with O" supported near the upper 
end of the string, Of supported near the lower end of 
the string, with the support of neither of those operators 
touching the boundary and with distance at least 2 be- 
tween the supports of those operators, and with O™ com- 
muting with the Hamiltonian. Again, the proof of this is 
very similar to lemma [2j So, physically in the toric code 
model one may say that these anyon excitations need to 
be connected by a string. 

Lemma[3]shows that these properties depend upon the 
existence of certain central elements in the interaction 
algebra. This is not surprising: we expect that if a theory 
has anyons then it is possible to make a measurement on 
B which detects the existence of an anyon inside B, and 
this measurement is precisely measuring the given central 
elements in the interaction algebra on B. 

Now we show the main result in this section on break- 
ability: 

Theorem 1. Define sets S and B as before. Define a 
set A which is an arbitrary set of sites. Suppose that 

A\I = XUY (44) 

for two disjoint sets X, Y such that for there is no pla- 
quette Z which has non-vanishing intersection with both 



X and Y. (A possible choice of such sets is shown in 
Fig.\@. 

Suppose Hb decomposes into Hb^e <£> T~Lb^i so that 
A E acts only on Hb->e cmd similarly A 1 acts only on 
T-Lb^tI- Then, any operator O supported on A that 
commutes with the term algebra is breakable on sets 
X U B,Y U B,I. 

Proof. Break An I into two disjoint sets, Ii,l2, such that 
h ny = h n A = 0. (For example, in Fig.SJ h is the set 
of sites in / above the middle of the heaviest solid line.) 
Such sets always exist as we could pick I\ = (A n I) \ Y 
and A2 = (A D I) \ I\ , though the choice in Fig. |4] is a 
different choice of Ji, I-i- Let A\ = XUli and A2 = VU^- 
Then we can write 

° = E°M 2 , (45) 

a 

with Oa 1 , Oa 2 supported on A\ , A2, respectively, where 
we can choose that 

a^f3 tr ((O^O^) =0 (46) 

and similarly and 

a^f3 -)■ tr((OS 2 ) f Ol) = 0. (47) 

This decomposition (|45[) can be found by a singular value 
decomposition. 

Note that 0% fulfills the the conditions of lemma [5] 
To show this, we need to show that the commutator 
[0% i: Qz} = for Z such that Z n E ^ 9. This com- 
mutator is vanishing if Z PI A\ = as the supports are 
disjoint. Suppose Z C\ A x ^0. Then, since Z H / = 0, 
ZnX ^ 0. Thus, ZHY = 0. So, [Qz,0 A2 ] = for all (3 
for such Z. So, since [O, Qz] = 0, usingthe orthogonality 
Eq. (|47f it follows that all the O a commute with those 

Qz also. Similarly, the 0^ 2 fulfill the conditions of the 
lemma also, with the set X replaced with Y. 

So, we can apply the lemma to both Al and 0^ 2 to 
get 

O = ^ OxubO B \j(a 1 \x)0 < b\j(a 2 \y)Oy\jb^ (48) 

a 

where 0^ uB and Oy uB both commute with the term 
algebra. Note that the index a in Eq. (|4"8")) is not the 
same as the index a in Eq. (|45|) ; we used Eq. ([45]) to 
write O as a sum of products of operators supported on 
A\ and A2 and then we used lemma[3]to write O a as sum 
of products of operators onlUB and B U (Ai \ X) and 
to similarly decompose 0% 2 as another sum of products 
of operators on Y U B and B U (A2 \ Y), thus writing 
O as a sum of products of operators on the four sets 
XUB,BU(A 1 \X),BU(A 2 \Y),YU B, as in Eq. 01 
Let Of = B u(a 1 \x)Obu(a 2 \y); so that 

= Y.°xubO i O YuB . (49) 

a 



We have thus written O as a sum of products of operators 
with the desired support, and Oxub and Or us both 
commute with the term algebra. So, if we could show 
that Of commutes with the term algebra, then we would 
have shown breakability. 

In fact, we cannot show that Of commutes with the 
term algebra in general, but by applying a singular value 
decomposition again we can fix this problem. Write the 
above equation as 

= J2°xuYu B Of, (50) 

a 

where 0^ uyuB = 0%ub°yub- b Y performing a sin- 
gular value decomposition as before, we can write this 
as 

= Y, 6 xuYu B Of, (51) 

a 

where now the operators XuYu b are orthogonal with 
respect to the Hilbert-Schmidt inner product and where 
Ox u yub i s som c linear combination of the operators 
XuYu b as follows: 

Oxuyub = ^ VctpO XyJ y uB , (52) 
P 

for some matrix V and similarly the operators Of are 
linear combinations of the operators Of . Now, using the 
orthogonality of OxuyuB' ^ follows that the operators 
Of commute with the term algebra. Plugging Eq. (|52[) 
into Eq. (|51[) . we find that 

= Y,V a pO xuB 0^ UB Of, (53) 

a/3 

so we have succeed in finding the desired decomposition 
of operator O; one may combine the two indices a, (3 into 
a single index to write this in the form Eq. <jTUJ). □ 

V. SOLVABILITY AND TRIVIAL STATES 
WITHOUT CENTRAL ELEMENTS 

Based on the discussion in section IIV1 where the ab- 
sence of central elements in a certain interaction algebra 
implies breakability, one might expect that the absence 
of these elements will also us to show that 2DCLH, re- 
stricted to instances without certain central elements, is 
in NP. We now show this (see below for a description of 
which interaction algebras we are considering). In the 
next section we consider instances with these central el- 
ements. In this section, where these central elements are 
absent, we are able to prove more than just that the 
problem is in NP: we also prove trivality (under a circuit 
definition of topological order using ancillas as discussed 
below) of a zero energy ground state in theorem [21 This 
triviality then shows that if a state is nontrivial under a 
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FIG. 5: Illustration of holes. The sites of the lattice are on 
the intersections of the solid horizontal and vertical lines. The 
heavier dashed lines represent the closed lines Bi. There are 
6 lines, corresponds to i = 1...6. Each such dashed lines has 
been drawn around a 4-by-4 square, so it contains 12 sites. 
For each line Bi, the 4 sites inside the square represent li. 

circuit definition (for example, a toric code ground state 
is nontrivial under this definition, even allowing the use 
of ancillas), then such central elements must be present. 

See Fig. [5j We have drawn several closed lines Bi , each 
line being the boundary of some set Si. 

Definition 8. A set of holes is a set of sets of sites, 
such that each set Si contains some set of four sites in a 
plaquette and such that Si DSj = for i ^ j . For each i, 
let Bi be the boundary of Si as defined previously and li 
be the interior. Let Ei be the set of sites not in Si . Let 
A\ be the interaction algebra of Si on Bi and let A\ be 
the interaction algebra on Bi of the algebra generated by 
the Qz for Z <f_ Si. 

Let be the Hilbcrt space on Bi. The algebras Af 
and A\ commute. So, if at least one of these algebras 
does not contain central elements, then we can decom- 
pose 'H.B i = HBi^Ei ® Hsi^ii so that Af acts only on 
T-Lsi^Ei and similarly A\ acts only on H-Si^u- 

Lemma 4. Assume we have a set of holes Si. Sup- 
pose for all i it is possible to decompose H-Bi as H-Bi = 
Hsi^Ei ®T~i-Bi-tii so that Af acts only on H-Bi^Ei and 
similarly A\ acts only on T-Lsi^u- Let He be the sum of 
Qz over all Z such that Z is not a subset of Si for any i. 
Let Hi be the sum of Qz over all Z such that Z is a sub- 
set of some Si . Then, H has a zero energy ground state 
if and only if both He and Hj have zero energy ground 
states. 

Proof. The "only if direction of the implication is ob- 
vious: since H = He + Hj and both He and Hi are 
positive semi-definite, H > He and H > Hi. 

The "if direction of the implication follows from the 
assumption that the T-Lst decomposes in the given way; 
this implies that the Hilbert space of the whole system 
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decomposes into H = He <£> Hi where He acts only on 
He an d Hj acts only on Hi so that the spectrum of H 
is the sum of the spectra of He and Hi. Here we have 
defined 

H! = (£){h u ®Hb^i)- (54) 

i 

□ 

We now show that given this decomposition of HBi, 
the problem of whether H has a zero energy ground state 
is in NP, under some assumptions on the sets Si. The 
Hamiltonian Hj is a sum over i of Hamiltonians Hi = 

SzcBiU/i Qz- 

Definition 9. We say that a set of holes forms a square 
lattice of size 0(1) if the following hold. Let the centers 
of the holes Bi form a square lattice with distance Thole 
between the centers, as shown in Fig. [3| and let each hole 
have diameter dhole- Let both dhole and rhole be 0(1). 

In Fig. [51 r hole = 6. 

The Hamiltonian Hj is a sum over i of Hamiltonians 
Hi = ^ZzaBui Qz- If the holes form a square lattice 
of size O(l), then each Hi acts on a Hilbert space with 
dimension and hence its spectrum can be deter- 

mined in time poly(D) by matrix diagonalization, and so 
whether or not Hj has a zero energy state can be deter- 
mind in time Npoly(D). Alsi, in this case the Hamilto- 
nian He can be coarse-grained into a Hamiltonian which 
is a sum of commuting two-body terms, with the di- 
mension of the Hilbert space in the coarse-grained prob- 
lem being D°^\ To do this coarse-graining, draw lines 
connecting the center of each hole to its four neighbors 
(above, below, left, and right of the given hole); these 
lines then divide the lattice of sites into different sets (de- 
form the lines slightly so that no site is exactly on a given 
line) ; then, coarse-grain by combining all sites in one set 
into a single supersite, and the resulting Hamiltonian is 
a sum of two-body terms (here we use the assumption 
that every hole contains at least one plaquette). Seei2 
for an illustration of a similar coarse-graining. Thus, the 
problem of whether He has a zero energy ground state 
is in NP, showing that the problem of whether H has a 
zero energy ground state is in NP. 

Note that both He and Hi have trivial ground states 
(trivial in the sense that it can be constructed by a local 
quantum circuit of depth and range O(l)); this result for 
He follows from the fact that He is a sum of commut- 
ing two-body terms on the supersites resulting from the 
coarse-graining. We claim that 

Theorem 2. Assume that HBi decomposes as above for 
a set of holes that form a square lattice of size O(l). 
Then, H has a trivial ground state, where we allow the 
use of ancillas in the construction of this ground state. 
That is on each site we are allowed to tensor in additional 
ancilla degrees of freedom, then initialize the system to a 
product state, then apply the quantum circuit, so that the 



resulting state should be a ground state of ' H®I , meaning 
H acting on the first copy tensored with the identity on 
the second copy. 

Proof. Define two copies of the system. On the first copy, 
called the real copy, construct a trivial state that is a 
ground state of He using a quantum circuit. Such a 
state follows from the coarse-graining above. On the 
second copy, called the ancilla copy, construct a trivial 
state that is a ground state of Hi. Then, apply another 
round of the unitary quantum circuit to swap the state 

on ^iiW-Bi-tii ® Hi^j between the two copies; this can 
be accomplished with a quantum circuit of depth 1 and 
range dhole] this circuit is a product of unitaries on the 
different B l U I t . □ 

It seems possible that this construction can be modified 
to avoid the use of ancillas, but we do not consider this 
here. 

We now describe an alternate way of showing that the 
problem of determining whether or not He has a zero 
energy ground state is in NP. This approach is based on 
the idea of propagation, as in Eq. (|T5|) . For the rest of the 
section, when we refer to the "term algebra" , we mean the 
algebra of terms of He ■ The set of sites in the column C 
which are not in any interior Ii can be written as a union 
of several intervals of sites , which we write X\ . q , X2 , c j • • • s 
with no interactions coupling sites in different intervals. 
Let the sets X^cs^G.c, ••■ be sets of sites in column C 
between the holes as shown in the figure. The minimal 
projectors which are central elements of the interaction 
algebra of the term algebra on C can be written as 

PZ 1 '"*- = P?-" P?-" .... (55) 

wth P'xfL supported on X a ,c- So, if there exists a zero 
energy ground state of He, then for each column of holes, 
we pick a column C that intersects the middle of that 
column of holes, and we can pick a set of a a ,c for every 
a so that TJ a c ^atf IS consistent with He- Note that 
ria ^atf IS a matrix product operator with dimension 
D°^\ So, by lemma |T]) the set of these matrix product 
operators is a witness for the existence of a zero energy 
state. 

VI. SOLVABILITY IN THE GENERAL CASE 

We now consider the case in which the decomposition 
used above does not hold. We actually begin with an- 
other special case: a case where the state may be topo- 
logically nontrivial but where certain density matrices 
are full rank and prove the existence of a witness in this 
case. We then discuss the general case. In a sense, one 
may regard the three cases we consider (the one in the 
previous section and the two in this section) as follows: 
the first problem considered is a problem which is es- 
sentially classical (no topological order) but where there 
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FIG. 6: One particular column C intersecting a column of 
holes. The solid lines represent the intervals X a ,c- In this 
case, this is column C — 4 and there are 4 intervals X a ,c for 
a= 1,...,4. 



such that the projector P^ 3 onto H.23 is in B^; indeed 
P23 is the operator that maps to the identity operator 
in the given full matrix operator under the isomorphism. 
The space T-L 2 3 is annihilated by all algebras B a . Up to 
a change of basis on ^ e isomorphism from the full 

matrix algebra to B23 is given by mapping a matrix in 
the matrix algebra to one or more copies of that matrix 
and the dimension of H.23 * s an integer multiple of the 
dimension of the matrices. 

Let the boundary algebra be the direct sum of alge- 
bras B23 over all a such that P^ does not annihilate the 
ground state subpace of Q2X3 ■ 

Let \i/jq) be a complete basis for the ground states of 
Q2X3 on sites 2,X, 3. Then, since the boundary alge- 
bra commutes with Q2X3, for every operator O in the 
boundary algebra, we have 



(58) 



may be "glassiness" ; the second problem considers a case 
where there may be topological order but where there is 
a certain assumption given below that one may interpret 
as having no glassiness; finally, we consider the case that 
combines both. 

A. Bound on Matrix Product Operator Dimension 
in 4-site Mode: Special Case 

In this subsection we consider a system of 4 sites, called 
1,2,3, X. Consider a Hamiltonian 

H = Qux + Q2X3, (56) 

where the operators Q12X and Q2X3 are positive semi- 
definite and Hermitian but need not be projectors. Let 
[Q\2XiQ2X3\ = 0. Let all sites have finite Hilbert space 
dimension, and let the Hilbert space dimension on site 
2 be D and let the Hilbert space dimension of site 3 be 
D3. Note: this quantity D3 will not appear in many of 
our bounds later, but we define this quantity because we 
will use it in intermediate steps of the calculation. 
We define 

Definition 10. Consider the algebra of operators on 
%2 <8>%3 generated by the full algebra of operators on site 
2 and by the interaction algebra of Q2X3 on site 3. Con- 
sider the subalgebra of this algebra consisting of operators 
which commute with the interaction algebra of Q2X3 on 
sites 2,3. 

Note that since this subalgebra is finite dimensional, it 
is a direct sum of subalgebras, each of which is isomorphic 
to a full matrix algebra. We call these subalgebras $23- 
The Hilbert space H.2 ® W3 can be written as a sum of 
Hilbert spaces, 

(0^2 3 )©^23 (57) 

a 



for some matrix X . We will often refer later to the matrix 
X corresponding to an operator O, by which we mean the 
matrix X that appears in Eq. Q58j) for the given operator 
O. Let X a be the set of such matrices X that correspond 
to operators O in B23', note that this set forms an alge- 
bra, since if operators O, O' G B^ 3 have corresponding 
matrices X, X', then the operator OO' is in £>2 3 also and 
has corresponding matrix XX'. Eq. (|58|) gives a homo- 
morphism from $23 to X a . 

Lemma 5. Suppose P^ does not annihilate the space of 
ground states of Q 2x3- Then, there is an isomorphism 
from B% 3 to X a . 

Further, let the effective boundary algebra be the 

direct sum of X a over a such that P^ does not annihilate 
the ground state subspace of Q2X3- Then, there is an 
isomorphism from the boundary algebra to the effective 
boundary algebra. 

Proof. Suppose a is such that P^ docs not annihilate 
the space of ground states of Q2X3- The algebra B% 3 
is isomorphic to a full matrix algebra. Suppose a given 
matrix M in the full matrix algebra maps to a given O 
in #23 under this isomorphism. Up to a change of basis 
which commutes with Q2X3, O must be some number 
of copies of M acting on the ground state subspace and 
some other number of copies of M acting on the rest 
of the space, padded with zeroes. However, since P^ 
docs not annihilate the space of ground states, there must 
be at least one copy of M acting on the ground state 
subspace. Thus, the matrix X that corresponds to O 
from Eq. (fB"8"|) must be, up to change of basis, one or 
more copies of M padded with additional zeroes. So, 
this gives the desired isomorphism from $23 to X a . 

We can similarly show the existance of an isomorphism 
from the boundary algebra, which is the direct sum of 
such B a , to the effective boundary algebra, which is the 
direct sum of such X a . □ 
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We will bound the dimension of the boundary alge- 
bra under certain assumptions. The first result, lemma[S] 
which is given in this subsection, will not be used later in 
the paper. However, the assumptions in this lemma are 
applicable to many physical models as discussed later, 
and so this lemma is worth mentioning. Also, this con- 
struction motivates the constructions of the next two sub- 
sections. The second result, lemma |H] which is given in 
the next subsection, is what we will use later. 

Lemma 6. Assume that H has a unique zero energy 
ground state ip° . Further, if the algebra on 2X generated 
by the interaction algebras oj Q\ 2 x an d Q2X3 on 2X has 
central elements P a projecting onto different subspaces 
ofHzx, then assume that there is some given a, which 
we write as a = 0, such that every state "0 that is an 
eigenstate of Q2X3 with zero eigenvalue also obeys P°ip = 
tp. The subspace H^X decomposes into two subspaces, 
H2X-ti <8>%2X-s-37 with Q12X acting only on %2X~>i and 
Q2X3 acting only on H.2X^>3- We require that Q2X3 has 
a unique zero energy state on / H 2 x^3®T~L 3 . Let K denote 
the rank of the reduced density matrix on TL 3 of this zero 
energy state. 

Let Dba be the dimension of the boundary algebra con- 
sidered as a vector space. For such a Qi 2 x,Q2X3, 

Dba < D 2 D 3 /K. (59) 

In particular, if the reduced density matrix is full rank 
so D 3 = K then Dba < D 2 and if D 3 /K = 0(1) then 
Dba < 0(1)D 2 . 

Proof. As discussed above, the subspace H% x decom- 
poses into two subspaces, H.2X^\ ^^2x^3, with Q2X3 
acting only on 0.2x^,3- Any ground state of Hamilto- 
nian Q2X3 is a pure state on %2Y->3 ® ^3- Let this 
pure state be <p. Introduce a basis |/i)2X->i f° r states on 
W^x-ti- S° the ground states |^°) above can be written 
as \n) 2 x^i ® <j>. 

We bound the dimension of the effective boundary al- 
gebra. Then, since there is an isomorphism from the 
boundary algebra (which is a direct sum of full matrix 
algebras) to the effective boundary algebra, this bounds 
Dba- 

Consider the vector space of operators acting on H2 ® 
H 3 . This space has dimension D 2 D 2 . Every operator 
in the boundary algebra is in this vector space. In fact, 
given any operator O in the boundary algebra, and any 
operator S acting on T~L 3 , the operator SO is in the space 
of operators acting on H2 ®T-L%. For such an operator O, 
we have that 

50|<> = sYt^M) (6°) 

= (j^»2X-n) ® (V). 

The above equation gives a map from the space of ground 
states of Q2X3 to the space spanned by states \y)ix^t\ 
tensored with arbitrary states on 7^2X^3 ® ^3 ( n °te that 



acting with an appropriate matrix S on <p can give an 
arbitrary state on %2Y->3 ® So, for each opera- 

tor SO, we have some corresponding linear map between 
these two spaces of different dimension. We will show 
that the dimension of the space of these maps is at least 
D 3 KDba- However, this implies that the space of op- 
erators on Ti.2 Cg> T-L 3 has dimension at least D 3 KDba- 
Hence, Dba < D 2 D 3 /K as claimed. 

To lower bound the dimension of the space of these 
maps, we define the Hilbcrt-Schmidt inner product of 
two different maps, one given by matrices O (and corre- 
sponding X from Eq. (|58[) ) and S and the other given by 
O' (and corresponding X') and 5". This inner product is 

YsWP^&S'O'^) (61) 

f 

= tr(X* X')(ct>\S f S'\<f>). 

Given an orthonormal set of matrices X using the inner 
product ti(X^X') and another orthonormal set of matri- 
ces S using the inner product (<j)\S^ S'\(/)), the product of 
these two sets gives us an orthonormal set of linear maps 
using the inner product (|6ip . Since the reduced density 
matrix of <j) has rank K on 3, we can construct D 3 K or- 
thonormal vectors with respect to that inner product. So, 
we can construct D 3 KDba different linear maps which 
are orthonormal using inner product (|6ip . □ 

Corollary 1. A similar bound on Dba holds also under 
the following slightly weaker assumptions. Suppose that 
the interaction algebra of Q2X3 on 3 is not the algebra of 
all operators but rather a subalgebra. Suppose the Hilbert 
space 14.3 can be decomposed into a product of two Hilbert 
spaces %\®1~L 2 , such that the interaction algebra of Qix3 
on site 3 only acts on the first space %\ in the product. 
So, Q2X3 can be written as some operator Q2X3 1 which 
is supported on 2,X and %\. We then define a new 4- 
site problem, replacing site 3 by the Hilbert space H\ and 
replacing Q2X3 by Q2X3 1 ■ Let %\ have dimension D 3 i . 
For this new A-site problem, if it turns out that the con- 
ditions of lemma\(j\ hold, with the state <j> having rank K 
on %\ then the bound Dba < D 2 D 3 i / K holds. 

While lemma [5] requires several assumptions about the 
Hamiltonian, these assumptions are in fact satisfied by 
many physical models after a certain grouping of sites. 
We now explain this grouping (which will be used in the 
proof of the main theorem later) , and use this grouping to 
show that there is a classical witness for the existence of 
a zero energy ground state for a certain class of 2DCLH 
problems (we explain which problems these arc after we 
explain the grouping). 

For each plaquette Z , define v(Z) to be the vertical 
coordinate of the bottom of the plaquette. Note that 
1 < v(Z) < L — l. 

Lemma 7. Consider an instance of 2DCLH. For each 
column C , for every i with 1 < i < L — 2, define the 
following A-site problem, which we call the A-site problem 
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centered at i. Group all sites in column C with vertical 
coordinate less than i into site 1. Group all sites in col- 
umn C with vertical coordinate greater than i into site 
3. Let site 2 be the site in column C with vertical coor- 
dinate equal to i. Group all sites in columns C — 1 and 
C + 1 into site X. Define Q\2X to be the sum of Qz 
over Z such that Z contains only sites which have been 
grouped into sites 1,2, X. Define Q2X3 to be a weighted 
sum of Qz over Z such that Z contains only sites which 
are grouped into sites 2, X, 3; these are the Z such that 
v(Z) > i. The weighted sum is a sum over such Qz, 
multipling each Qz by some generic positive scalar; the 
reason for this multiplication is so that every operator 
which commutes with Q2X3 will commute with each such 
Z separately; note that the zero energy states of Q2X3 co- 
incide with the zero energy states of an unweighted sum. 

Suppose that for each C and i, the resulting A-site prob- 
lem obeys the conditions of lemma\$for some given K, D3 
(different values of K,D 3 may occur for different C,i). 
Let (D3/ K) max denote the maximum of D3/K over all 
C,i. 

Then, for each column C , every central element Oc of 
the interaction algebra of C — 1,(7, C + 1 on C can be 
written as a sum of two terms 

c = 0™ ps + 0h, (62) 

where O mps is a matrix product operator with dimen- 
sion bounded by D^{Dz/K) max and O^Pc-i.c+i = 
Pc-i,c+iOq = 0. Further, if the given central element 
of the interaction algebra is a projector, then both O mps 
and O are projectors. 

Proof. Take the central element Oc- For each vertical 
coordinate i, let 1_ba(*) be the unit of the boundary al- 
gebra for the 4-site problem centered at i. Then, for 
any 3, 2 < j < L - 2, let 0™ ps (j) = IbaU^baU ~ 

1)...1ba(2)O c 1ba(2)...1baO' - Let 0™ ps = 

0™ ps (L-2). 

To show that this operator 0™ ps has the desired prop- 
erties, consider the 4-site problem centered at j. The 
operator 0™ ps (j) commutes with Q2X3 in the given 4- 
site problem, because each of the operators lBA(k) for 
k < j commute with the given Q2X3- For k = j, this 
follows by definition of the boundary algebra. For k < j, 
it follows because IbaQ^) commutes with the Q2X3 in 
the 4-site problem centered at k and Q2X3 in the 4-site 
problem centered at j appears with a generic coefficient 
in the weighted sum defining Q2X3 in the 4-site problem 
centered at k. 

Hence, 0™ ps (j) can be written as a sum of products 
of operators: 

OcH^EOt-O?, (63) 

where Of is supported on site 1 after grouping (sup- 
ported on sites in column C with vertical coordinate 
less than j) and O® is supported on sites 2, 3 after 
grouping and Of commutes with Q2X3 in the 4-site 



problem centered at j. Further, because 0™ ps (j) = 
1ba(j)0™ ps (i)lsA(j)) the operator Of is in the bound- 
ary algebra. Since we have bounded the dimension of the 
boundary algebra by D 2 (D 3 /K) max , we have bounded 
the number of terms in this sum. 

Write 0™ ps (j) as a matrix product operator; we can- 
not bound the bond dimension on all of the bonds, but 
we can bound the bond dimension on bonds connecting 
site i to i + 1 for i < j. This follows inductively since 
if 0™ ps (j — 1) obeys this bound, then each Of is a ma- 
trix product operator with bounded bond dimension. So, 
0™ ps is a matrix product operator with bond dimension 
at most D 2 (D 3 /K) max . 

Further, 0™ ps (j) — 0™ ps (j — 1) vanishes on the ground 
state of Q2X3, so OqPc-i,c+i = as claimed. 

Finally, since 1ba(j) is a projector and 1ba(J) com- 
mutes with 0™ ps (j — l) (this follows because 1ba{j) is 
in the center of the algebra of operators in sites 2, 3 that 
commute with Q2X3), it follows that O mps (j) is a pro- 
jector. □ 

The conditions of the above lemma with (D 3 / K) max = 
0(1) actually hold for many physical models, such as 
the toric code and Levin- Wen models with appropriate 
boundary conditions, such as the case in which the pla- 
nar model is given by "squashing" a cylinder into a plane 
(although of course for that model we know the exis- 
tence of a zero energy state by other methods). Assume 
a zero energy state exists a model with such a bound 
on (D3/ K) max . Here is a first attempt, which unfortu- 
nately fails, at constructing a witness to the existence of 
this state. Given that a zero energy state exists, there 
exists a sequence of projectors Oc, one projector for each 
column O, which are minimal central elements of the in- 
teraction algebra such that tr(FJ c Oc Yl z Pz) > 0. For 
each projector Oc, we write Oc = Oq PS + Oc- It seems 
at first as if the operators 0™ ps meet the conditions of the 
operators Pc in lcmma[TJ Unfortunately, even though we 
know that Oc commutes with Pc-i,c, we don't know 
that this holds for 0™ ps , preventing the use of lemma [1] 

There are several ways around this problem. We will 
give one approach here. The following lemma is almost 
the same as lemma [TJ we omit the proof since it is the 
same except for a slightly different grouping, grouping 
only column C — 1 in X, not C — 1 and C + 1 and ex- 
cept for the use of corollary [T] instead of lemma [S] Note 
that we could have chosen to use corollary [T] instead of 
lemma [6] in lemma [7] we chose not to do it there since 
the conditionsof lcmma|6]apply to many interesting prob- 
lems. However, for the next lemma, for many problems 
lemma [6] would not give a useful bound but corollary 
Q] still does; for example, for the toric code Hamiltonian, 
D 3 i j K would be of order unity but D3 / K would not have 
any useful bound with the grouping of the next lemma. 

Lemma 8. Consider an instance of 2DCLH. For each 
column C , for every i with 1 < i < L — 2, define the 
following 4-site problem, which we call the A-site problem 
centered at i. Group all sites in column C with vertical 
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coordinate less than i into site 1. Group all sites in col- 
umn C with vertical coordinate greater than i into site 
3. Let site 2 be the site in column C with vertical coor- 
dinate equal to i. Group all sites in columns C — 1 into 
site X . Define Q12X to be the sum of Qz over Z such 
that Z contains only sites which have been grouped into 
sites 1,2, X. Define Q2X3 to be a weighted sum of Qz 
over Z such that Z contains only sites which are grouped 
into sites 2,X, 3. The weighted sum is a sum over such 
Qz, multipling each Qz by some generic positive scalar; 
the reason for this multiplication is so that every operator 
which commutes with Q2X3 will commute with each such 
Z separately; note that the zero energy states of Q2X3 co- 
incide with the zero energy states of an unweighted sum. 

Suppose that for each C and i, the resulting A- site prob- 
lem obeys the conditions of corollary [7] for some given 
K,D 3 i (different values of K, D 3 i may occur for different 
C,i). Let (-D31 /K) max denote the maximum of D 3 i/K 
over all C, i. 

Then, for each column C , every central element Oc of 
the interaction algebra of C — 1, C on C can be written 
as a sum of two terms 

c = 0™ ps +0h, (64) 

where O mps is a matrix product operator with dimen- 
sion bounded by D 2 (D 3 i /K) max and OqPc-x,C+x = 
Pc—x,C+xOq = 0. Further, if the given central element 
of the interaction algebra is a projector, then both O mps 
and O are projectors. 



Then, it follows that 

Theorem 3. Consider an instance of 2DCLH which has 
a zero energy ground state. Suppose the conditions of 
lemma\^hold, for (D\/ K) rnax < poly(N). Then, there 
exists a witness to the existence of that ground state which 
can be efficiently verified. 

Proof. Given that a zero energy state exists, there exists 
a sequence of projectors Oc, one projector for each col- 
umn C, which are minimal central elements of the inter- 
action algebra such that tr(]T c Oc Y\ z Pz) > 0- For each 
projector Oc, write Oc = 0™ ps + Oq, with 0™ ps ,Op 
as in lemma |5J The projector 0™ ps commutes with 
Pc-i.c (this follows from the fact that Oc commutes 
with Pc~x,c an d that OqPc-x,c = for this approach). 
Also, 0™ ps is a mask as follows from the fact that Oc 
is a mask and that O^Pc-i.c = 0- So, by lemma [TJ the 
set of projectors Oc form a witness to the existence of a 
zero energy state. □ 

Note that the conditions of this theorem hold for many 
interesting Hamiltonians such as Levin- Wen and toric 
code (though of course for those models we know the 
existence of zero energy ground states by other means). 



B. Bound on Matrix Product Operator Dimension 
in 4-site Mode: More General Case 

We now return to the 4-site model, but consider dif- 
ferent assumptions. 

Lemma 9. Consider the A-site model. Suppose that there 
is no projector P in the boundary algebra such that the 
partial trace ^(P^xs-f ) * s not full rank on 3 , where P2X3 
is the projector onto the zero energy subspace of Q2X3 ■ 

For such a Qx2X, Q2X3, for each a such that P a does 
not annihilate the ground state subspace of Q2X3, the di- 
mension of the boundary algebra, considered as a vector 
space, is bounded by Dba < 

Proof. Wc bound the dimension of the effective boundary 
algebra, as in the proof of lemma [5] 

Consider the vector space of operators acting on H2 <S> 
H3. This space has dimension D 2 D\. Every operator 
in the boundary algebra is in this vector space. In fact, 
given any operator O in the boundary algebra, and any 
operator S acting on H3 , the operator SO is in the space 
of operators acting on H2 <8> "H^. For such an operator O, 
we have that 

S0\^)=SJ2X^\%)- (65) 

V 

The algebra of matrices X is a direct sum of one or more 
algebras X a which are isomorphic to full matrix algebras 
Ai a . Consider the algebra M a for some given index a 
and let it have dimension D a . For each such algebera 
Ai a , consider it as a vector space and consider the sub- 
space of the vector space corresponding to matrices which 
are non-zero only in the first row (of course, this choice 
is somewhat arbitrary, as we could make an arbitrary 
change of basis). Let V a be the vector space of matrices 
in X a that maps to this given vector space under the iso- 
morphism. Choose a basis for V a consisting of matrices 
that map to a matrix with a 1 in the z-th column of the 
first row and a zero everywhere else, for each 1 < % < D a . 
Call the matrices in this basis X(i,a). 

We define the Hilbert-Schmidt inner product of two 
different maps, one given by matrices O (and correspond- 
ing X) and S and the other given by O' (and correspond- 
ing X') and S'. This inner product is 

£(<|oWo W) (66) 

Note that if X = X(i,a) and X' = X(j,/3), then this 
inner product is zero unless i = j and a = (3. Further, 
for each i,a, we can set X = X' = X(i,a) and use 
Eq. ([66]) to define an inner product for the S, with D| 
orthonormal basis vectors. This is the point at which we 
use the assumption that the reduced density matrix on 
site 3 of P2X3P f° r every projector P in the boundary 
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algebra is full rank, since the inner product for S, is 
tr 3 (SS^ p^), for some density matrix p 3 which is a reduced 
density matrix of such a P^xz-P- 

So, this gives us ^ Q D a D 2 different choices of X, S 
which are orthonormal using inner product (|66[) . So, 
E Q D a < D 2 . This implies that J2 a D l < P> 4 - □ 

We have the following generalization of this lemma, 
which is slightly different to corollary [TJ using instead a 
direct sum decomposition: 

Lemma 10. A similar bound on Dba holds also under 
the following slightly weaker assumptions. Suppose that 
the interaction algebra oj Q2X3 on 3 commutes with some 
projector P3 on H3. Let the dimension of the range of 
P3 be D 3 i . Suppose for that all ground states of Q2X3 
are in the range 0/P3. Construct the A-site problem as 
above. For this new A-site problem, if it turns out that 
the conditions of lemma hold, with the state <f> having 
rank K on W3 then the bound Dba < D 2 D 3 i / K holds. 

Proof. Since the interaction algebra commutes with P3, 
we can write it as a direct sum of an algebra which does 
not annihilate P3 and an algebra that does annihilate P3 . 
Every operator in the second algebra annihilates every 
vector in the ground state of Q2X3- The dimension of 
the algebra that does not annihilate P3 is at most D 2 t . 
We then use the same proof as in lemma [5] □ 

C. Reduced Dimensionality in General Case 

The above lemma leads to the following result in the 
case of a general instance of 2DCLH. Consider a general 
instance for which a zero energy state exists. Then, there 
exists a sequence of projectors Oc , one projector for each 
column C, which are minimal central elements of the in- 
teraction algebra such that tr(J7 c Oc Y\ z Pz) > 0. Sup- 
pose that these projectors cannot be written as low bond 
dimension matrix product operators. Using lemma |H] or 
[TU1 a result analogous to lcmma[5]can be proven, showing 
the existence of a decomposition Oc = Oq + 0™ ps , as- 
suming that the conditions of lemma |H] or [TU] are met for 
every grouping for some D 3 i/K which is at most poly(N) 
(we abbreviate this by simply saying that the conditions 
of the lemma "are met" below, rather than every time 
referring to the ratio D 3 i/K). So, if Oc is not a low 
bond dimension matrix product operator, either the con- 
ditions of lemma [TU] are not met for some grouping, or 
1 ba is not equal to the identity operator for some group- 
ing. We claim that in either case, we can construct a 
projector II supported on the set of sites in column C 
with vertical coordinate greater than or equal to j, for 
some j such that II commutes with all the Qz except 
possibly for those Qz with v{Z) = j, and such that II is 
consistent with the Hamiltonian, and such that II is not 
full rank. 

To prove this claim, suppose first that 1 ba is not the 
identity operator for some grouping. Then, simply take 



II = 1b a- Instead, suppose that the conditions offTOlarc 
not met for some grouping centered at i. Then, there 
is some projector R in the boundary algebra such that 
its reduced density matrix on sites i + 1, i + 2, ... is not 
full rank. Suppose that there is a ground state such 
that (ip, Rip) 7^ 0. In this case, we will take II to be the 
projector onto the range of this reduced density matrix. 
We now show that IT has the desired properties. Let 
X be the plaqucttc in columns C — 1,C with v(Z) = i 
and let Y be the plaquette in columns C, C + 1 with 
v(Z) = i. Then, QxQyRiP has non-zero norm, as may be 
verified by noting that \QxQyR^\ > ', Q 'xQy Rip)\ = 
\{ip,Rip)\. Further, QxQyR^ is a ground state, as R 
commutes with all Qz except for Z = X and Z = Y, 
and HQxQyR^' = QxQyRiP, so II is consistent with 
the Hamiltonian. The commutation requirements on IT 
hold also. 

Suppose instead that no such ground state ijj exists. In 
this case, we take II = 1 — R. Thus, in either case, such 
a projector exists. 

Using these projectors IT leads to the following idea 
for a proof in the general case. We treat even and odd 
columns differently, constructing the groupings and pro- 
jectors Oc only for even C. We define a new problem, 
H' = ^2 Z Q Z , wn ere in a deviation from our previous 
notation each projector Q' z is supported on the set of all 
sites in the same column as a site in Z with a vertical 
coordinate greater than or equal to v(Z) — 1 (that is, the 
sites in Z as well as all sites "below" them). Given two 
projectors A,B, define A n B to project onto the inter- 
section of the range of A and B, and define A U B to 
equal 1 — (1 — A) n (1 — B). So, we follow the following 
iterative procedure. Start by defining Q' z to be cup over 
all Qy such that Y is in the same two columns as Z and 
v(Y) > v(Z). Then, support of the Q' z is indeed as we 
claimed. Then, if the operators Oc are not low bond di- 
mension matrix product operators, construct a projector 
IT as discussed above, and replace Q' z = Q' z U (1 — IT) 
for all Z such that Z intersects column C and such that 
the support of IT does not include any sites above Z. 
This new problem still has a zero energy ground state, 
because IT was consistent with H. We repeat this con- 
struction again and again for different groupings until no 
more such IT can be found. 

Now, if this new Hamiltonian H' were still defined by a 
set of commuting projectors, this would allow us to find 
a witness for H, because using results above we could 
show that central elements of H' are low bond dimension 
matrix product operators. (The fact that the support of 
H has been increased will not affect our use of lemma [UJ 
since wc really only used properties of Q2X3 there; the 
4-site problem for H' now has terms Q12X3 and Q2X3 
but we only use properties of the latter term). However, 
the new Hamiltonian may not be a set of commuting 
projectors: the projectors Q' x and Q' Y need not commute 
if v (X) = v(Y) and X and Y intersect on an even column. 
This may occur because n may not commute with Qz- 
We will explain elsewhere how to deal with this situation. 
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VII. DISCUSSION 

We have proven the existence of classical witness to 
the existence of a zero energy ground state for certain 
instances of 2DCLH; a general proof will be given else- 
where. Other results in sections IlIIIIVIVl help character- 
ize possible central elements in the interaction algebra in 
these systems. 

One interesting feature is that we naturally write cer- 
tain central elements in the interaction algebra as ma- 
trix product operators. We have also shown how to do 
this specifically for certain central elements in Levin- Wen 
models. This representation as a matrix product opera- 
tor may have useful practical applications in attempts to 
engineer a system that realizes a Levin- Wen model, per- 
haps using Josephson junctions or cold atoms. At first 
glance, realizing a Levin- Wen model seems completely 
impractical, since one needs to realize twelve spin oper- 
ators. However, the ability to represent these operators 
at matrix product operators may help lead to a practical 
realization by interacting given degrees of freedom with 
a low dimensional auxiliary system; this leads to an in- 
teresting question of physically implementing a matrix 
product operator with bounded bond dimension. We do 
not consider this question further. 

The bound on the dimension of the boundary algebra 
may be interesting for use elsewhere. This bound can be 



interpreted as a bound on the number of particle types 
possible in a topological quantum field theory described 
by a lattice model. After all, to create a pair of anyons 
in such theories, one acts with a stringlike operator that 
commutes with the Hamiltonian everywhere except at 
its ends. This stringlike operator does not annihilate the 
ground state. Assuming this operator is in the algebra 
generated by the interaction algebra and by the algebra 
of operators at its endpoints, this stringlikc operator is 
then in the boundary algebra. 

We relate the existence of certain central elements to 
the circuit definition of topological order in section [Vj 
showing triviality of the ground state of H in section [V] 
in the absence of certain central elements. 

Of course, a natural question is what happens in three 
dimensions. Is the problem of whether a sum of com- 
muting projectors with local interactions on a three di- 
mensional cubic lattice has a zero energy ground state a 
problem that is in NP? The discovery of codes like Haah's 
cod o 19 ' 20 in three dimension casts some doubt on whether 
this can be true. These codes show that there might be 
some very "wild" possibilities in three dimensions which 
make the decision problem much harder. 
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